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SUMMARY 

Statistics  of  the  tra jector ie-  generated  by  a  second-order  target 
model  are  presented,  and  using  these  statistics  the  distribution  of  the 
change  in  target  heading  over  any  time  period  is  derived.  This  distribu¬ 
tion  is  the  basis  of  a  new  technique  for  choosing  the  model  manoeuvre 
parameter  which  is  suitable  for  a  given  class  of  targets.  Design  curves 
for  this  technique  are  supplied. 
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1  INTRODUCTION 

An  essential  component  of  most  target  tracking  filters  is  a  mathematical  model 
describing  the  kinematics  of  the  target.  The  target  model  represents  the  level  of 
manoeuvre  which  the  target  of  interest  is  expected  to  exhibit,  and  is  used  in  the  track¬ 
ing  filter  for  smoothing  the  estimate  of  the  target  state  {eg  position  and  velocity). 

This  Report  describes  the  properties  of  a  simple  linear  second-order  target  model,  which 
is  widely  used  in  tracking  filters*.  The  chief  objectives  of  this  Report  are  to  present 
the  statistics  of  the  target  trajectories  described  by  the  target  model  and  using  these 
statistics  to  derive  a  new  method  for  selecting  the  manoeuvre  parameter  of  the  model 
appropriate  to  a  given  class  of  targets.  The  trajectory  statistics  are  also  used  to 
obtain  a  search  region,  such  that  there  is  a  given  probability  of  a  target  which  obeys 
the  second-order  model  being  within  the  region  at  some  time  after  observing  the  target 
state  {ie  assessing  how  quickly  targetting  information  becomes  stale). 

The  continuous  version  of  second-order  target  model  is  derived  under  the  assumption 
that  the  acceleration  of  the  target  in  each  Cartesian  co-ordinate  is  an  independent  Gaussian 
white  noise  process.  The  discrete  version  of  the  model  is  similar  except  that  the  target* 
acceleration  is  taken  to  be  a  Gaussian  white  sequence,  the  acceleration  being  assumed 
constant  over  each  time  step.  The  model  has  the  advantage  of  simplicity  while  providing, 
for  many  applications,  quite  an  adequate  representation  of  target  trajectories  (in  two 
or  three  dimensions).  The  trajectory  described  by  the  model  is  a  variation  about  a  con¬ 
stant  velocity  course,  whose  magnitude  and  direction  are  defined  by  initial  conditions. 

The  deviation  from  this  mean  course  is  controlled  by  a  single  parameter,  the  power 
spectral  density  of  the  Gaussian  driving  noise  (this  parameter  is  taken  to  be  the  same 
for  each  Cartesian  co-ordinate  in  this  Report).  If  this  parameter  is  zero,  then  the 
target  trajectory  is  a  constant  velocity  course. 

The  main  contribution  of  this  Report  is  the  derivation  of  the  probability  density 
function  (pdf)  of  the  change  in  heading  angle  ( ie  the  change  in  direction  of  the  velocity 
vector)  in  a  given  period,  for  a  target  whose  motion  is  governed  by  the  second-order 
model.  This  pdf  is  a  function  of  the  model  manoeuvre  parameter,  the  initial  speed  of  the 
target  and  the  time  period  over  which  the  heading  change  is  considered,  and  it  provides 
a  useful  measure  of  the  degree  of  target  manoeuvrability  which  is  being  modelled.  Hence 
this  distribution  of  heading  change  forms  the  basis  of  a  new  method  for  selecting  the 
manoeuvre  parameter  appropriate  to  a  given  class  of  targets  (section  4).  In  section  3 
the  problem  of  defining  a  search  region  for  a  target  obeying  the  second-order  model  and 
given  prior  imperfect  targetting  information  is  discussed.  Assuming  the  manoeuvre  para¬ 
meter  of  the  target  is  known,  it  is  shown  that  a  suitable  search  area  for  the  two- 
dimensional  case  is  an  ellipse  (or  an  ellipsoid  for  three  dimensions)  centred  on  the 
expected  position  of  the  target.  The  search  region  and  the  distribution  of  heading  angle 
change  are  both  derived  from  the  conditional  pdf  of  target  state,  which  follows  directly 
from  well-known  general  results  for  linear  systems  driven  by  Gaussian  noise.  These 
general  results  are  collected  in  Appendix  A,  which  also  discusses  the  relationship  between 
linear  systems  and  the  Markov  property. 
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Save  for  the  derivation  of  the  model  in  section  2,  the  chief  emphasis  of  this  Report 
is  on  the  discrete  version  of  the  target  model,  since  nowadays  tracking  filters  are 
usually  implemented  digitally.  However  at  the  end  of  sections  3  and  4  and  Appendix  A, 
the  chief  results  are  also  presented  for  the  continuous  version  of  the  target  model. 

2  DERIVATION  OF  THE  TARGET  MODEL 

The  continuous  version  of  the  target  model  is  based  on  the  premise  that  the  accel¬ 
eration  of  the  target  is  a  zero  mean,  Gaussian,  white  noise  process.  Hence  each  Cartesian 
co-ordinate  x  may  be  described  by* 

x  -  w’  ,  (1) 

where  w*  is  a  Gaussian  noise  process  with 

E[w*]  -  0  1 

and  >  ■  (2) 

E[w'(t  +  r)w'(t)]  «  q’S(t)  j 

where  E  is  the  expectation  operator,  6(t)  is  the  Dirac  delta  function  and  q'  is  the 
power  spectral  density  (constant)  of  the  acceleration  noise.  The  white  noise  process  w’ 
may  be  viewed  as  the  limit  of  a  bandlimited  Gaussian  noise  process  with  a  constant  power 
density  spectrum  q'  up  to  the  cut-off  frequency  .  The  autocorrelation  function  of 
such  a  bandlimited  process  is 

sin  ui  t 

E[v'(t  ♦  T)v’(t)]  -  q' - 


and  as  ®  ,  so  the  noise  becomes  white, 

sin  u  t 

■  7TT^~  " 

(see  Papoulis*^,  p  98),  from  which  equation  (2)  follows. 

The  general  state  space  representation  of  a  continuous  linear  system  driven  by 
Gaussian  white  noise  is  as  follows: 

x  -  Fx  ♦  Gw'  ,  (3) 

where  x  is  the  state  vector,  F  and  G  are  matrices  which  may  be  time-dependent  and 

*  A  rigorous  statement  requires  the  use  of  Ito  stochastic  calculus,  in  which  case 
equation  (I)  would  be  written 

dx  »  du  , 

4 

where  du  is  a  Wiener  process  with  zero  mean  and  variance  q'dt  (see  Sage  and  Melsa 
or  Jazwinsky^  for  further  details). 
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w'  is  a  vector  of  white  Gaussian  noise  (sometimes  called  process  noise), 
it  is  assumed  that 

E[w*]  *  0 

and  ,  p 

E|w'(t  +  T)w'^(t))  -  Q'  6  (t) 


In  this  Report 

W 


where  Q'  is  the  power  spectral  density  matrix  of  the  process  noise. 
The  target  model  can  be  put  into  this  general  form  by  setting 


Q’  ■  q’  . 


(5) 


The  state  transition  matrix  (see  Ref  5)  of  this  system  is 

/'  T 

$' (t  ♦  T,t)  -  I 

\0  I 

which  depends  only  on  the  time  interval  i  .  It  can  be  seen  that  the  target  velocity 

.  2 

x  is  the  integral  of  w'  and  is  therefore  a  Wiener  process  (see  Papoulis  ). 

For  the  discrete  version  of  the  target  model,  the  target  acceleration  is  assumed 
to  be  a  white  Gaussian  sequence,  the  acceleration  being  constant  over  a  time  step. 

Hence,  integrating  over  a  time  step  At  ,  for  each  Cartesian  co-ordinate, 

x  j(k  ♦  l)At)  -  x(kAt)  ♦  Atx(kAt)  ♦  jAt^w^  ] 
and  v  ,  (6) 

x  j(k  ♦  l)At)  -  x(kAt)  ♦  Atw^  I 

where  (w^)  is  a  white  Gaussian  sequence  of  variance  q  ,  ie 

EK'  ' 


0 


6 


and 


E(V'j1 


q5 


kj  * 


The  general  state  space  representation  of  a  discrete  linear  system  driven  by  a  Gaussian 
white  noise  sequence  is 

*k+i  ’  +  r?k  •  (7) 


where  is  c^e  state  vector  at  time  kAt  , 

is  the  state  transition  matrix  from  k  to  k  +  1  , 
T  is  a  matrix 

and  { }  is  a  random  Gaussian  sequence  with 


E  [-k‘ 


and 


:M] 


Q5kj  * 


where  Q  is  the  covariance  of  the  sequence  {w^  }  . 

The  discrete  target  model  given  by  (6)  can  be  written  in  this  general  form  as  follows: 


1 1  a  A 

(Mzn\ 

5k+l 

\o  1  j 

K  +  1 

\  At  / 

(8) 


x(kAt) 


f  x(kAt)\ 
^x(kAt)/ 


and  Q  *  q  . 

For  this  model,  the  target  velocity  x(kAt)  is  a  random  walk  ^see  Jenkins  and  Watts'*)  . 
If  the  value  of  q  is  chosen  so  that 

q  ■  fr  ’ 

it  can  be  shown  that  as  At  -►  0  and  the  samples  became  dense,  the  discrete  target 
model  (8)  tends  to  the  continuous  model  (5).  To  see  this,  ewrite  equation  (8)  as 

“<wv  ■  (.  .K’l.K 

and  let  At  -»  0,  k  -  «  and  j  -*  -  such  Chat  kAt  ■  C  and  jAt  *  t  .  In  this  case  it 

4 

can  be  shown  that  (see  Sage  and  Melsa  ), 
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and 


E  [w.  w .  ] 

*  j 


q6. 


kj 


At  *kj  "q'6(C 


T) 


which  is  the  continuous  model. 

3  THE  PROBABILITY  DENSITY  FUNCTION  OF  THE  TARGET  STATE 
3. 1  The  discrete  case 

3.1.1  Propagation  of  the  mean  and  covariance  of  the  state  vector 

The  state  space  representation  of  the  target  model  for  discrete  time  steps  is 
given  by  equation  (8),  and  from  Appendix  A,  section  A. 1  it  can  be  seen  that  the  sequence 
of  state  vectors  resulting  from  this  model  is  a  homogeneous  Gauss-Markov  sequence.  Suppose 
it  is  known  that  the  pdf  of  5^  is  Gaussian  with  mean  and  covariance  ,  but  that 

is  unknown.  From  equations  (A-9)  and  (A-10)  of  Appendix  A,  section  A. 1 ,  it  can  be 
shown  that,  after  an  interval  of  lit  ,  the  pdf  of  x^  conditioned  on  this  information, 
le  P^Sk+ll!Sk,}*k)  ’  *S  c’auss*an  “ith  mean 


Vt  "  *% 


\\ 4  (iAt)ik] 

A, 


(9) 


and  covariance 


“k 


+4  *  ♦  q(At)2 

i-o 


[t  ('  +2U-  i  -  I))]  f('*2 


(l  -  j  -  1) 


At 
\  2 


(l  ♦  2(1  -  j  -  I)) 


f  Ctnl  I >k  +  21At<m12>k  4  (lAt)  ♦  *At(m,,),. 


+  ^Y2~  (4!-2  '  1)1 


22  k  v  1 2yk  ‘“tv“22'k 

+  q(^t)3i2 
2 


3.2 


(“,2>k  4  iit(n22>k  4  ^T1 


(B22>k  4  q<lt)  1 


.  (10) 


8 


"  ^  * 

(m22)k/ 

of  the  covariance  matrix  is  the  variance  of  the  target  position  x  , 
is  the  variance  of  the  target  velocity  x  and  m^  is  the  correlation  between 
the  target’s  position  and  velocity.  Thus  equations  (9)  and  (10)  describe  the  propaga¬ 
tion  of  the  mean  and  the  covariance  of  the  target  state  vector  given  that  the  state 
vector  has  a  mean  of  and  a  covariance  matrix  .  Also  with  l  *=  1  , 

equations  (9)  and  (10)  are  implicit  in  the  a  -  S  tracking  filter  (in  which  each  target 
co-ordinate  is  filtered  separately),  equation  (10)  giving  the  covariance  of  the  target 
state  immediately  before  a  measurement  is  taken. 


/<n  n)k 

where 

\(Dl2)k 
The  element 


If  the  position  xk  and  velocity  of  the  target  are  known  perfectly,  then 

the  conditional  pdf  2^)  is  Gaussian  with  mean 


(\  *  UAt)\\ 

V  *k  I 


(I!) 


and  covariance 


q(4trn 


1  (At)2(4Z2  -  I) 
12 


Ut 

2 


l  it 
2 


(12) 


These  relations  may  be  derived  from  (9)  and  (10)  by  letting  -*■  0  so  that  x^  , 

or  directly  from  (A-5)  and  (A-6)  of  Appendix  A,  section  A.l.  If  the  initial  target  state 
xQ  is  known  perfectly,  then  the  expected  path  of  the  target  is  a  constant  velocity 
trajectory  (see  (11))  and  the  probable  deviation  from  this  path  is  given  by  the  con¬ 
ditional  covariance  (12).  It  can  be  seen  that  the  conditional  covariance  ( 1  2)  is  only  depen¬ 
dent  on  the  interval  l  ,  and  the  conditional  variance  of  position  increases  approximately 
with  (ZAt)  ,  while  that  of  velocity  increases  with  (ZAt)  . 

3.1.2  Likelihood  ellipses  of  the  Gaussian  pdf  of  the  state  vector 

The  propagation  of  the  pdf  of  x  may  be  visualized  in  -he  state  space  (x,x)  by 
the  evolution  of  cor  .ours  of  constant  probability  density  or  likelihood.  For  the  Gaussian 
pdf  of  x  ,  these  contours  are  ellipses  and  are  known  as  likelihood  ellipses.  The  centres 
of  these  evolving  ellipses  lie  on  the  mean  state  space  trajectory  which  is  a  constant 
velocity  path  represented  by  a  straight  line  parallel  to  the  x-axis  in  (x,x)  space. 

Suppose  chat  the  initial  state  vector  Xq  is  Gaussian  distributed  with  covariance  Mq  . 

If  the  position  is  initially  uncorrelated  with  the  velocity,  ie  j 2 ^ 0  "  then  the 
likelihood  ellipses  corresponding  to  x^  will  have  their  principal  axes  aligned  with  the 
x  and  S  axes  of  the  state  space.  Now  suppose  that  there  is  no  system  driving  noise 
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(q  -  0)  so  that  the  target  trajectories  will  have  constant  velocities  equal  to  their 
initial  velocities.  In  this  case  the  likelihood  ellipses  at  a  later  time  £At  may  be 
obtained  via  a  linear  shear  transform  of  the  initial  ellipse: 
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Thus  the  likelihood  ellipses  will  become  slewed  with  increasing  time  (see  Fig  1)  as 
targets  with  initial  velocities  greater  than  the  mean  will  travel  further  than  the  mean 
target,  while  targets  with  initial  velocities  less  than  the  mean  will  travel  less  far. 

In  other  words,  the  correlation  between  the  position  and  velocity  increases  with  time 
and  this  effect  is  described  mathematically  by  the  first  term  in  the  covariance  propaga¬ 
tion  equation  (10),  ie 

l  *T 

m£  -  <P  M0* 


If  the  system  driving  noise  is  non-zero  (ie  q  >  0) ,  then  the  target  trajectories  are  no 
longer  constant  velocity  paths  and  this  introduces  further  uncertainty  into  the  target 
state  with  increasing  time,  so  increasing  the  area  of  the  ellipses  (see  Fig  2).  This  is 
represented  by  the  second  term  of  equation  (10),  which  is  superposed  with  the  first 
effect . 

If  the  initial  target  state  vector  is  known  perfectly  (ie  MQ  *  0) ,  then  as  already 
discussed,  the  propagation  ot  the  mean  and  covariance  of  the  Gaussian  pdf  of  the  state 
vector  is  given  by  equations  (11)  and  (12).  For  this  case,  in  Figs  3  to  7 ,  the  evolution 
of  the  likelihood  ellipses  has  been  shown  for  several  different  values  of  the  acceleration 
noise  variance  q  (the  manoeuvre  parameter).  In  these  diagrams,  the  velocity  axis  has 
been  normalized  to  the  mean  velocity  v^  (which  is  also  the  initial  target  velocity)  and 
the  distance  axis  has  been  noramlized  to  VgAt  ,  the  mean  distance  travelled  in  one  time 
step.  The  initial  state  is  set  to 


-52  -12 

and  q  is  varied  between  10  (v^/At)  and  10  (v^/At)  ,  the  value  increasing  by 

a  factor  of  10  in  each  Figure.  In  each  case  likelihood  ellipses  of  the  state  vector  are 
shown  for  l  -  10,  30,  60  and  100.  The  ellipses  have  been  chosen  so  that  there  is  a  95Z 
probability  of  the  state  falling  inside  the  ellipse.  Also  shown  in  each  Figure  is 
a  realization  of  the  state  space  trajectory  generated  by  driving  equation  (8)  with 
a  random  Gaussian  sequence  of  variance  q  ,  and  the  four  points  marked  by  squares  on  the 
trajectory  are  the  states  corresponding  to  the  four  likelihood  ellipses.  In  almost  every 
case  these  points  are  within  the  95%  likelihood  ellipses. 
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3.1.3  Search  regions  for  the  second-order  target  code! 

It  is  often  necessary  to  define  a  search  region  for  a  target  based  on  prior 
imperfect  targetting  information.  If  the  target  obeys  the  second-order  model  given  by 
equation  (8),  and  the  uncertainty  in  the  targetting  information  for  each  Cartesian 
co-ordinate  is  independent  and  can  be  expressed  in  terms  of  the  covariance  matrix 
(which  may  be  different  for  each  co-ordinate),  then  the  pdf  of  the  target  position,  at 
a  time  £At  after  the  observation,  is  Gaussian  with  each  co-ordinate  independent  and 
having  a  mean  given  by  (9)  and  a  variance  (m^)^  given  by  (10).  If  the  search  region 
is  defined  as  that  region  for  which  there  is  a  given  probability  of  the  target  being 
within  the  region,  and  for  which  the  value  of  the  position  pdf  at  every  point  within  the 
region  exceeds  the  value  of  the  pdf  outside  the  region,  then  in  two  dimensions  the  search 
area  is  bounded  by  a  likelihood  ellipse  and  in  three  dimensions  it  is  bounded  by  an 
ellipsoid  (see  Bryson  and  Ho^) ,  with  the  expected  target  position  at  the  centre  of  the 
region.  In  two  dimensions,  the  search  area  is  defined  by 


.  f  (x>y)  <  i  . 

mn  °n  2  2 

x  y 


where  x  ■=  xQ  +  ZAtxQ  (9), 


mn  ’  +  2iit(°I2 v)n  *  WAt)2(ta22v)n  +  (4*2  _  >>l  fre" 


‘VO 


x  0  12 


(10), 


is  the  observation  of  the  x  co-ordinate  position  and  velocity, 


/<*ll,>n  <m12Jr 


V®i2x)0  <"22  xy 


and  y  and  n>]  are  similar. 


is  the  covariance  of  the  error  in  this  observation, 


Likewise,  in  three  dimensions  the  search  volume  is  defined  by 

(x  -  x)2  t  (y  -  ;>2  t  ( 2  -  z)2  -  f 3 (x ,y , z)  <  T3 


1  1 


Since  f-(x,y)  is  the  sum  of  the  squares  of  two  independent  Gaussian  variates,  it  has 

2  7  .  2 

a  x  distribution  with  2  degrees  of  freedom,  and  similarly  f^Cx.y.z)  has  a  x  dis¬ 
tribution  with  3  degrees  of  freedom.  Hence  if  the  required  probability  of  the  target 
being  within  the  search  region  is  P  ,  then  the  thresholds  T2  and  may  be  found 

from  tables  of  the  distribution.  For  example,  for  P  -  0.95,  “  6.00  and 

T3  -  7.81  . 
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1  ] 

3 . 2  The  continuous  case 

In  this  brief  section,  the  equivalent  continuous  versions  of  the  major  results  cf 
section  3.1  are  given.  These  results  may  be  derived  by  setting  q  =  q'/lt  and  letting 
At  -*■  0  and  l  ■*  »  so  that  ZAt  -*•  t  . 

The  stochastic  process  x(t)  generated  by  the  continuous  target  model  (5)  is 
a  homogeneous  Gauss-Markov  process.  If  at  time  t  ,  the  pdf  of  x(t)  is  Gaussian  with 
mean  x(t)  and  covariance  M(t)  ,  then  the  pdf  cf  x(t  +  t)  conditioned  on  this  infor¬ 
mation,  ie  p|x(t  +  x)  |*(t)  ,M(t))  ,  is  Gaussian  with  near. 


x(t)  +  ix(t)\ 

k  *(t)  j 


and  covariance 


4  THE  DISTRIBUTION  OF  CHANGE  IN  TARGET  HEADING  AND  A  METHOD  FOR  SELECTING  THE 
MANOEUVRE  PARAMETER 

4 . I  Discrete  case 

In  this  section,  the  pdf  of  the  change  in  heading,  over  a  given  period  of  time,  of 
a  target  obeying  the  second-order  model  is  given.  The  details  of  the  derivation  of  this 
pdf  are  given  in  Appendix  B.  The  distribution  of  heading  change  is  the  basis  of  a  method 


12 


I* 


for  choosing  a  suitable  manoeuvre  parameter  q  for  a  given  class  of  targets.  The 
information  required  for  choosing  q  is  an  estimate  of  the  mean  speed  of  the 
targets  and  the  likely  maximum  change  of  target  heading  in  some  time  period. 


First  consider  a  target  moving  in  two  dimensions  with  each  Cartesian  co-ordinate  of 
the  target  governed  by  the  second-order  target  model  (equation  (8))  with  manoeuvre  para¬ 
meter  q  .  The  change  in  the  heading  angle  of  the  target  cculd  be  positive  or  negative, 
but  since  it  is  the  magnitude  of  the  change  which  is  of  interest,  only  the  modulus  of 
the  change  in  heading  angle  is  considered.  In  Appendix  B  the  pdf  of  the  change  in 
heading  angle  8  after  a  period  lAt  (where  0  takes  values  in  the  range  [O,-*))  has 
been  derived: 


p(9) 


03) 


where 


cos  0 


and  erf(.)  is  the  error  function  defined  by 


erf  (x)  -  —  [ 


The  distribution  function  of  8  defined  by 


F(0)  *  I  p(A)d> 


X-0 


is  the  probability  of  the  changa  in  heading  angle  being  less  than  3  ,  and  does  not  appear 
to  have  any  simple  analytical  representation.  F(8)  may  be  interpreted  as  the  probability 
that  the  target  vector  remains  within  a  sector  of  semi-angle  0  with  axis  of  symmetry 
along  the  original  target  heading  (see  Fig  8a).  It  can  be  seen  that  p(6)  and  F(0) 

depend  only  on  the  parameter  cl  and  are  independent  of  the  initial  target  heading. 

-2  -I  2 

Figs  9  and  10  show  p(0)  and  F(0)  for  ci  taking  the  values  10  ,10  ,1,  10,  10  and 

103.  From  equations  (13),  as  ;2  ■*  “  so  p(0)  -*■  1/*  and  F  ( 0 )  0/*  (also  see  Figs  9 

and  10),  indicating  that  for  large  values  of  C  or  after  a  large  number  of  time  steps 
the  distribution  of  9  becomes  uniform  (ie  there  is  no  preferential  target  heading)  - 
as  would  be  expected. 

To  make  use  of  the  distribution  of  0  for  the  selection  of  a  suitable  value  of  q 
for  modelling  a  class  of  targets,  contours  of  constant  value  of  F ( 0 )  have  been  plotted 
in  Fig  1 1  as  a  function  of  9  and  the  distribution  parameter  Cl  .  Contours  are  shown 
for  F (9)  taking  the  values  of  0.3,  0.6,  0.9,  0.95  and  0.99,  for  the  distribution  para¬ 
meter  in  the  range  [10~3,  I02J  .  If,  for  the  class  of  targets  of  interest,  the  value 


TR  85071 


13 


of  9  fcr  which  F(9)  was  equal  to  one  of  the  contour  values  were  known,  then  the 
required  value  of  the  distribution  parameter  (which  determines  q)  could  be 

directly  read  from  Fig  II.  However  such  precise  information  is  unlikely  to  be  available, 
although  there  is  good  chance  of  the  tracking  filter  designer  having  some  estimate  of  the 
likely  maximum  change  of  heading  within  some  time  period.  If  this  likely  maximum  change 
of  heading  is  formally  defined  as  the  value  of  0  for  which  there  is  a  95Z  probability 
of  the  target’s  change  of  heading  being  less  than  0  ,  then  this  value  may  be  associated 
with  the  F(9)  ■  0.95  contour  of  Fig  11.  Hence  if,  for  the  class  of  targets  of  interest, 
the  maximum  likely  heading  change  in  a  time  £At  is  approximately  0'  ,  then  a  suitable 
value  of  q  for  modelling  these  targets  is  given  by  (see  definition  of  $  with 
equation  (13)) 


where  (C£)’  is  read  from  Fig  11  as  the  value  corresponding  to  F(9*)  in  the  region  of 
0.95.  For  example  consider  a  class  of  targets  with  a  mean  speed  of  about  15  m/s  with 
an  expected  maximum  heading  change  of  */ 4  in  13  seconds  (corresponding  to  a  turning 
circle  of  radius  250  metres).  From  Fig  11,  F  (7t/4  )  -  0.95  corresponds  to  a  value  of 


C —  0.15  ,  which  indicates  that  a  suitable  value  of  q  would  be  about  (for 
At  -  1  second). 


To  illustrate  the  type  of  trajectories  generated  by  different  manoeuvre  parameters, 

2 

Figs  15  to  19  show  realizations  of  target  trajectories  for  q  *  C(Vg/At)  where  ; 
takes  the  values  10  10  \  10  \  10  ^  and  10  *.  Each  Figure  shows  six  realizations  of 

a  100  time  step  trajectory  (solid  line),  as  well  as  the  expected  trajectory  (dotted  line) 
and  the  95Z  likelihood  ellipse  for  the  last  point  in  the  trajectory  (given  the  exact 

initial  conditions).  When  considering  the  change  in  heading  angle  over  the  hundred  tine 

-3  -2  -1 

steps  (£  »  100),  the  parameter  takes  the  values  10  ,10  ,10  ,1  and  10  for 

Figs  15  to  19  respectively.  For  the  admittedly  limited  sample,  the  observed  change  of 
heading  is  in  agreement  with  the  theoretical  distribution  (Fig  11).  For  example  in 
Fig  17  for  which  Z.I  ■  0.1  ,  only  one  of  the  9ix  trajectory  realisations  shows  a  heading 
change  as  much  as  0.5  radian  (F(0.5)  -  0.9)  and  for  the  remaining  examples  the  change  is 
less  than  0.25  (F(0.25)  -  0.6). 


The  pdf  of  the  change  in  heading  angle  after  i  time  steps  has  also  been  derived 
for  a  target  moving  in  three  dimensions  (see  Appendix  B) : 


2a  f 

'  ^  |5 


(1  +  200  ♦  erf(O)}  sin  8 


where  £ 

and  C 


cos  9 


► 


14 

distribution  function  may  be  simply  expressed 
_ l_ 

2ci  2 

— j t-  (l  ♦  erf  (?))  cos  8  .  (15) 

In  three  dimensions,  F (8)  may  be  interpreted  as  the  probability  that  the  target  velocity 
vector  remains  inside  a  cone  (of  infinite  extent)  of  semi-angle  9  ,  with  apex  at  the 
origin  and  axis  of  symmetry  along  the  original  target  heading  (see  Fig  8b).  As  for  the 
two-dimensional  case,  p(0)  and  F ( 9 )  are  independent  of  the  initial  target  heading 
and  depend  only  on  the  distribution  parameter  £2  (p(9)  and  F(e)  are  shown  in 
Figs  12  and  13  for  5 1  taking  values  of  10  ^ ,  10  1,  10,  10^  and  10^).  As 

equations  (14)  and  (15)  indicate  that 


Also,  unlike  the  two-dimensional  case,  the 
(see  Appendix  B) : 


F(0) 


*  ('  *  "'fef  ■ 


and 


P  W 


sin  6 
2 


06) 


Since  the  proportion  of  the  total  volume  of  space  enclosed  by  the  infinite  cone  of  semi- 

2 

angle  9  increases  as  sin  (9/2)  ,  the  distribution  function  (16)  corresponds  to 
a  velocity  vector  pdf  which  has  spherical  symmetry  about  the  origin  (see  Appendix  B  for 
a  more  rigorous  argument).  Hence,  as  would  be  expected,  as  £2  -*•  •  there  is  no  prefer¬ 
ential  target  heading.  Fig  14  shows  contours  of  F(9)  taking  the  values  of  0.3,  0.6, 

—3  2 

0.9,  0.95  and  0.99  for  £2  in  the  range  (10  ,10  1  .  This  Figure  can  be  used  for 

choosing  i  suitable  value  of  q  for  modelling  targets  moving  in  three  dimensions  in  the 
same  way  as  for  the  two-dimensional  case. 

The  following  points  concerning  this  technique  for  selecting  q  should  be 
considered : 

(i)  The  technique  can  only  be  sensibly  employed  for  expected  maximum  heading  angle 

changes  in  the  range  (0.1,2.8]  radians,  which  corresponds  approximately  to  £2  varying 

-3  ... 

between  10  and  1.  Outside  this  range,  a  small  change  in  maximum  heading  angle  results 

in  very  large  changes  in  the  parameter  £2  and  since  expected  maximum  heading  change  is 

unlikely  to  be  known  very  precisely  this  could  lead  to  large  errors  in  the  estimate  of  q  . 

(ii)  For  the  class  of  targets  of  interest,  an  estimate  of  the  maximum  likely  heading 
change  in  1  time  steps  will  probably  be  available  for  several  different  values  of  2  , 
and  for  each  of  these  a  corresponding  value  for  q  can  be  obtained.  If  the  values  of  q 
do  not  differ  by  more  than  one  order  of  magnitude,  then  it  is  suggested  that  the  average 
value  of  q  should  be  taken.  However,  if  the  values  of  q  vary  over  several  orders  of 
magnitude,  the  second  order  model  may  be  an  inadequate  representation  of  the  class  of 
targets  (ie  the  random  walk  model  for  target  manoeuvre  is  incorrect)  and  it  may  be 
necessary  to  consider  an  alternative  model. 
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(iii)  In  selecting  the  value  of  q  on  the  basis  of  the  maximum  likely  change  in 
heading  (and  using  the  second-order  model),  there  is  an  implicit  assumption  that  the 
distribution  of  heading  angle  change  for  the  class  of  targets  of  interest  is  similar  to 
the  corresponding  F(0)  . 

(iv)  The  value  of  q  could  be  determined  via  some  other  statistic  of  the  change 
in  heading  angle  if  this  information  were  available  rather  than  an  estimate  of  the  likely 
maximum  change.  For  example,  for  certain  applications  the  mean  change  of  heading  angle 

9  after  some  period  might  be  available.  To  determine  q  from  this  information,  the 
relationship  between  0  and  would  have  to  be  obtained  via  numerical  integration 

using  equations  (13)  or  (14). 

(v)  This  technique  for  selecting  q  can  only  be  expected  to  give  an  order  of 
magnitude  estimate.  However,  this  is  probably  adequate  for  many  applications  since  the 
performance  of  tracking  filters  is  sensitive  to  the  order  of  magnitude  of  q  rather  than 
to  the  exact  value,  and  also  the  lack  of  information  (or  uncertainty)  in  the  likely 
manoeuvres  for  a  target  class  will  often  preclude  a  more  precise  determination  of  q  . 

If  further  information  on  the  target  trajectories  likely  to  be  encountered  is  available, 
it  may  be  possible  to  make  a  more  precise  determination  of  q  by  simulation. 

4.2  Continuous  case 

The  equivalent  manoeuvre  parameter  for  the  continuous  second-order  model  is  q’  , 
which  is  the  constant  power  spectral  density  of  the  acceleration  noise  (see  section  2). 
The  pdf  of  change  of  target  heading  in  x  seconds  for  the  continuous  case  is  also  given 
by  equations  (13)  and  (14)  for  the  two-dimensional  and  three-dimensional  cases  respect¬ 
ively,  if  the  distribution  parameter  is  replaced  by  nr  ,  where  n  is  given  by 

2.-2. 

n  *=  q'/Vg  and  has  units  of  s  /Hz  (if  the  power  spectral  density  q'  is  with  respect 
ro  Hz).  Hence  Figs  11  and  14  may  also  be  used  for  choosing  the  manoeuvre  parameter  for 
the  continuous  model.  If  the  value  of  (Ct)  appropriate  to  the  known  likely  maximum 
heading  change  in  time  x  is  read  from  Figs  11  or  14,  as  indicated  in  section  4.1,  then 
a  suitable  value  of  q*  for  the  continuous  model  is  given  by 


q* 


S£il  v2 

t  0 


5  CONCLUSIONS 

The  statistics  of  the  trajectories  generated  by  a  second-order  target  model  have 
been  presented.  The  mean  or  expected  target  trajectory  is  a  constant  velocity  path. 

The  uncertainty  (ie  the  variance)  in  the  future  position  of  the  target  about  the  expected 
path,  given  perfect  knowledge  of  the  current  position  and  velocity,  increases  approximately 
as  the  cube  of  time.  The  variance  of  the  target  velocity  increases  linearly  with  time. 

Both  the  position  and  velocity  variance  are  also  proportional  to  the  acceleration  noise 
variance  q  in  the  discrete  formulation,  or  the  constant  power  spectral  density  q’ 
in  the  continuous  case. 

Using  the  statistics  of  the  target  trajectories,  the  probability  density  function  of 
the  change  in  target  heading  over  a  given  time  period  has  been  derived.  This  distribution 


16 


forms  the  basis  of  a  new  technique  for  choosing  the  manoeuvre  parameter  q  which  is 
able  for  a  given  class  of  targets.  Design  curves  are  supplied  which  can  be  used  for 
selecting  q  ,  given  knowledge  of  the  likely  maximum  change  of  target  heading  in  some 
period  and  the  mean  target  speed. 


suit- 
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(A-9) 
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and 


Vi 


**V 


/'n'lrQrT>:'n'1 


(A-  1  0) 


The  above  equations  describe  the  propagation  of  the  mean  and  the  covariance  of  the 
sequence.  Note  that,  as  would  be  expected,  (from  (A-3) ,  (A-4) ,  (A-7)  ar.d  (A-8))  when 
*  0  ,  indicating  that  x^  *  , 

P(WVV  “  P<5k+t|5k)  • 


Finally  the  correlation  matrix  between  members  of  the  sequence  is  given  by 

Vi  k  -  E[(*k+Z  -  Vi><5fc  -  VT] 

=  ♦(k  +  Z.k)^ 

and 

\  k+1  -  vT(k  +  4-k>  • 

A. 2  Continuous  systems 

The  definitions  and  relations  for  continuous  systems  are  analogous  to  those  for 
discrete  systems.  A  continuous  random  (or  stochastic)  process  x(t)  is  Markov  if,  for 
any  t,  <  t2' <  . ..  <  tQ 

p(5(tn)l5(tn.,).S(tn_2)I...,*(t,))  -  p(5<tn>|x(tn_,))  . 

A  Markov  process  is  completely  determined  statistically  in  the  interval  (t^t^)  if  the  density 
functions  p(x(t j) |x(t2))  and  P(x(Cq))  are  known  for  all  t^  <  t^  in  the  interval 
(tp.t^).  A  Markov  process  is  homogeneous  if  pjx(t  ♦  T)|x(t)J  is  independent  of  t  . 

A  Gauss-Markov  process  is  a  Markov  process  for  which  p|x(t)  |x(r))  and  p^x(t)^ 
are  Gaussian  pdf.  The  state  vector  of  a  continuous  linear  system  driven  by  Gaussian 
white  noise  is  a  Gauss-Markov  process,  ie  if 

x  *  Fx  ♦  Gw'  ,  (A-ll) 

where  x  is  the  state  vector,  F  and  G  are  matrices  and  w'  is  a  vector  of  white 
Gaussian  noise  such  that 

E  (w'l  -  0 

and 

E[w'(t  ♦  *)v,T(t)]  =  Q'5(t)  , 

then  the  output  x(t)  of  this  system  is  a  Causs-Markov  process. 
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Expressions  which  are  analogous  to  the  discrete  case,  may  be  derived  for  the 
conditional  pdf  of  output  x(t)  of  continuous  system  (A-ll)  (see  Bryson  and  Ho^) : 

p(x(t  +  x)|x(t))  ■/[xlt  +  t)  ;  *  ,C  ’)  , 

u'  *  ♦  ' (t  ♦  T,t)x(t) 


where 

and 


C' 


j  <X> *  ( t  ♦  T,a)GQ'C^$,T(t  +  T,a)da 


(A- 12) 


(A- 13) 


where  V(tj,t.)  is  the  transition  matrix  of  (A-ll). 

If  $’(t  T,t)  is  only  dependent  on  t  ,  ie 

❖  *(t  ♦  x,t)  =  4>'(x)  for  all  t 


then  (A- l 2 )  and  (A- 13)  reduce  to 


k'  -  ♦'(i)x(t) 


and 


i 

C'  •  J  *'(x  -  a)GQ'GT{i'T(T  -  a)da  . 


(A-U) 


(A- 15) 


In  this  case  it  can  be  seen  that  p(x(t  +  t)|x(t)]  is  independent  of  t  and  so  the 
Gauss-Markov  process  is  homogeneous. 

If  the  pdf  of  x(t)  is  known  to  be 

p(x(t))  =.r(x(t);J(t),M(t)) 

then 

p(x(t  +  t)  |  x(t)  ,M(t))  -  „>^jx(t  ♦  t)  ;x(t  ♦  t),M(t  +  t))  , 

where  x(t  +  t)  -  »'(t  +  t,t)x(t)  (A-16) 


1.  ▼  \ 

M(t  +  T)  -  $'(t  ♦  t,t)M(t)4,T(t  ♦  t,t)  +  j  *'(t  +  T,o)GQ'GTi,T(t  +  t,n)do  (A-17) 


which  is  the  solution  of 


HM  T  T 

~  -  FM  ♦  MF  ♦  GQ’G‘ 

d  t 
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Finally  the  correlation  matrix  cf  the  process  is  given  by 

R(t  +  T,t)  =  O’ (t  +  T,t)M(t) 

and 

R(t,t  +  T)  =  M(t)»,T(t  +  T,t) 
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Appendix  B 

DERIVATION  OF  THE  DISTRIBUTION  OF  CHANGE  IN  TARGET  HEADING 
Two-dimensional  case 


In  this  Appendix,  the  pdf  of  the  change  in  heading,  over  a  given  period  of  time, 

for  a  target  obeying  the  second-order  model  (equation  (8)  of  section  2)  is  derived. 

Suppose  that  a  target  moving  in  the  x-y  plane  has  a  velocity  [*0]  at  some  instant. 

»  •  \  v  yo • 

The  velocity  at  a  later  time  has  a  pdf  which  nay  be  determined  from  the  results 

of  section  3.  If  each  co-ordinate  is  independent  and  each  has  the  sane  manoeuvre  para¬ 
meter  q  ,  then,  for  the  discrete  target  model,  equation  f ' 2)  of  section  3  indicates 
that  each  velocity  component  has  a  Gaussian  distribution  with  variance  q(lt )“I  and 
mean  x_  or  yn  after  a  time  interval  of  i^t  .  Hence 


p(*.y!vV 


_±_  «„[■-  Jjf,  -  v 


(y 


V 


(B-l) 


where  a  =  q(Ac)  Z  . 

The  change  of  heading  X  ,  is  the  angle  betueen  the  vectors  and  j  ^0^  .  Since 

only  the  absolute  value  of  the  change  of  heading  is  of  interest,  define  the  probability 
distribution  function  F(6)  as  the  probability  that  '  is  within  the  range  [-S.f]. 
Hence  F(9)  is  given  by 


'O)  ■  I!  P(i,yii0.y0)<lidy  , 


( B— 2 ) 


where  D  is  the  region  of  the  (x,y)  plane  for  which  Ac  [-9,3]  . 

To  write  the  integral  ( B— 2 )  in  a  form  more  amenable  to  analysis,  consider  a  change 

2.2.2 

of  variables  from  (x,y)  to  (r,A)  where  r  =  x  ♦  y  and  A  is  as  defined  above.  Ising 
the  cosine  rule  (see  Fig  Bla),  the  argument  of  the  exponential  in  (B-l)  may  be  written 


2  2  2  2 
(x  -  xQ)  +  (y  -  yQ)  =  r  +  r0  -  2rrQ  cos  A 

where  rQ  -  / i2Q  *  y2Q  . 

Hence  (B-2)  becomes 


F(e) 


2rr^  cos  \) 


rdrdA 


(3-3) 


Since  the  integrand  is  symmetric  with  respect  to  A  ,  changing  the  range  of  integration 
over  A  to  (0 , 9 I  and  doubling  the  value  of  the  right-hand  «ide  of  (B-3)  leaves  the  value 
of  F(9)  unchanged.  If  also  an  integration  by  parts  on  the  integral  over  r  is  per¬ 
formed,  then  (B-3)  becomes 
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o 

■i] 


I  L_r°/2a2  +  r° cos  x 


L ri5 1 


r  +  ro  •  2rro  cos  V  drfd;  • 


The  remaining  integral  over  r  can  be  evaluated  (see  Abramowitz  ar.d  Stegun  ,  7.4.2)  and 


2  2  . 
■ro/2a  ? 


i  +  /i£e^(I  +  erf(0)d*  , 


rQ  cos  7 

where  ^  =  - 

/2a 

Now  if  the  initial  speed  of  the  target  was  ,  so  that  vq  “  rQ  *  anc^  ^  c-  bitten 


where  £  is  a  non-dimensional  parameter,  then  (3-4)  cay  be  written 


F(6)  = 


f  r-  r2 

-j—  j  1  *  >^Ce'  (1  +  erf(0)d- 


where  £  =  —  . 

/ia 

Hence  the  distribution  function  of  heading  change  depends  only  or.  the  parameter  . 

Using  numerical  integration  F(9)  has  been  plotted  for  different  values  of  ;r.  in 
Fig  30.  The  pdf  of  heading  change  is  given  by 

p<e>  - 


and  so  directly  from  (B-5) 


where  £  = 


2;M  r2  1 

— - 1  I  *  /?-t  (1  ♦  erf  (())>  , 


p(0)  is  shown  in  Fig  9  for  different  values  of  the  parameter  Cl  . 

B. 2  Three-dimensional  case 

By  a  simple  extension  of  (B-l),  the  conditional  pdf  of  velocity  for  the  three- 
dimensional  case  is  given  by 


i.y.z'VW  ’  7r!r3  ex?f  73  [(i 

(2^)7  a  l  2a  *- 


The  probability  distribution  function  F(v)  for  this  . ..  : 

that  the  target  velocity  vector  remains  within  a  ccr.e  -- 
angle  S  ,  with  apex  at  the  origin  and  axis  of  symmetry  a  I 
F(9)  can  therefore  be  derived  by  integrating  (3-7)  over  t:.<  v  . 
facilitate  this  consider  a  spherical  co-ordinate  system  with 


from  the  vector  Yq  I  »  azimuthal  angle  measured  with  rcspc 


contains  the  vector  y^  and  radial  distance  r  given  by 


2  .2  .2  .2 

r  =  x  +  y  +  z 

From  the  cosine  rule,  the  argument  of  the  exponential  in  CB-7)  can  be  written  (see 
Fig  Bib) 

(x  -  x0)  ♦  (y  -  yQ)  +  (z  -  zQ)  =  r  ♦  rQ  -  -rr0  cos  -  , 

.  f.2  T  72 

whore  r0  =  JXq  *  Xq  +  zq  . 

Also  for  a  spherical  co-ordinate  system,  the  elemental  volume  is 

2 

r  sin  >  drd,dV 

Hence  the  distribution  function  of  heading  change  is  given  by 


- — j  I  I  exp j - -  |r2  +  r2  -  2rrQ  cos  r“  sin  a  d*d r  J  d;  .  (B-S) 

(2V’a  rf0  kiQ  l  2a'  1  d0 


The  integration  with  respect  to  >  can  be  performed  by  inspection,  and  with  evaluating 
the  integral  over  ;  ,  this  gives 


— 

r^Va  J.  r0 
r -0 


rrQ/a2  rr  coSeV\ 

:  -  e  dr 


(B-9)  can  be  written  as  the  sum  of  two  integrals  of  the  same  form  as  the  integral  over  r 
in  (B-3)  and  hence  can  be  integrated  in  the  same  manner  to  give 
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cos  -X 


f  (5)  =  ill  +  erf  — - 


e7  (I  +  erf(-))  cos  v  , 


If  the  initial  target  speed  was  ,  so  that 


it  q  is  written  as 


*  m  HTT 


F(0)  = 


i(i  •  *"fc 


2r  l  " 2 

- - ^  (1  +  erf  (O)  cos  0  , 


where  £  = 


Hence  the  distribution  function  only  depends  on  the  parameter  ; l  ,  and  in  Fig 
is  shown  for  different  values  of  ; 1 '  .  The  density  function  p(r)  may  be  obt: 
differentiating  (B-10)  to  give 


13  Fir) 


p(6)  = 


£  +  ~  e **  (1  +  2 •; “ )  ( 1  +  e r f  ( O )>  sin.  0 


where  K  = 


p(v)  is  shown  in  Fig  12  for  different  values  of 
As  -'2-  ■*  ®,  s  ■*  0  and 


F(9)  -  sin2(4) 


This  is  the  distribution  of  elevation  angle  which  results  from  any  pdf  with 
symmetry.  To  prove  this  consider  any  function  f  whose  integral  over  all 
unity  and  which  only  depends  an  r  .  Hence,  for  any  spherical  co-ordinate 


n  ® 

IS  I 

#«-0  x  =  0  r*0 


f(r)r  s in  >  drd>d  r  =  1 


which  implies 


j  f(r)r2 


:6 
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4  IT  ® 

/  d*  / 

i^=0  r=0 


£(r)r  dr 


-  (1  -  cos  0)2* —  =  sin  ( —  I 

4  it  \  21 

using  (B-12),  which  provides  Che  required  proof.  So  as  would  be  expected,  as  £2  -*►  ® 
there  is  no  preferential  heading. 

B. 3  Alternative  derivation 

An  alternative  derivation  of  the  distribution  of  change  in  heading  angle  using 
Student*s  non-central  t  distribution  has  been  suggested  by  David  Lloyd**  of  AW’2  Division, 
RAE  Farnborough.  If  t*  is  defined  by 


where  x 


<$  is  some  constant 

and  x,Xj > . . . ,x  are  independent  and  Gaussian  with  zero  mean  and  unit  variance, 
then  t'  has  the  pdf 


P ( t '  ;v,  <5) 


AJrU  Mr  2 


tfp) '  / "" '  i(-7r#rr) 


(which  follows  from  (9.6.7)  of  Fisz^).  (B-14)  is  an  expression  for  the  non-central 
t  distribution. 

First  consider  the  three-dimensional  case.  Rotate  the  Cartesian  co-ordinate  system 
about  the  origin  to  form  a  new  Cartesian  system  (x'jy'.z')  with  z'  aligned  along  the 
initial  target  heading  (the  orientation  of  x’  and  y*  is  arbitrary  from  symmetry 

2 

considerations) .  Then  xf ,  y'  and  z’  are  independent  and  Gaussian  with  variance  a  , 
and  x*  and  y'  have  zero  mean  while  z'  has  a  mean  Vq  (which  is  the  initial  target 
speed).  The  change  in  heading  9  is  the  angle  between  the  V  axis  and  the  vector 


/  .  ,  2  .,2 

v  x  ♦  y 

i"  ♦  (vQ/a) 


J.,,2  T  -,.2 

V  x  ♦  y 
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2? 


where  x",  y"  and  z"  are  independent  and  Gaussian  with  zero  nean  and  unit  variance 

g 

(see  Ruben  for  the  general  n  dimensional  case).  Comparing  this  expression  with 
(B— 13),  it  can  be  seen  that  /2  cot  0  has  a  non-central  t  distribution  with  v  *  2 
and  6  -  v^/a  .  Hence,  the  pdf  of  0  is  given  by 


J  t  * 

[>(••4 

-V-Zacote 


t’  =  /2  cot  6 


cosec  9 


Noting  that  t’  ■  /2  cot  0  implies  that 


sin^9  and 


and  putting  (B-14)  into  (B-15), 


oo 

/ 

/  2 

1 

/  u  exj>l 

" 2  r 

0  i 

\ 

2  follows 

using 

F(9) 

1 

•  1 

Likewise  for  the  two-dimensional  case,  the  change  in  heading  v  is  given  by 


where  x"  and  y”  are  Gaussian  with  zero  mean  and  unit  variance,  so  that  cot  6  has 
a  non-central  t  distribution  with  v  »  1  and  5  *  v^/a  •  Hence 


*-  Jt'-cot 


ao 

- 

2\2" 

f 

1 

2  v°  vo 

/  w  exp 

2 

w  -  2  —  w  cos  0  ♦  — =■  1 
a  2  I 

0 

3 1 J 

from  which  (B-3)  of  section  B.l  follows  using  (B-16). 
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